High-gradient operators in perturbed Wess-Zumino-Witten field theories in two 

dimensions 
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Many classes of non-linear sigma models (NLcrMs) are known to contain composite operators with 
an arbitrary number 2s of derivatives ("high-gradient operators") which appear to become strongly 
relevant within renormalization group (RG) calculations at one (or fixed higher) loop order, when 
the number 2s of derivatives becomes large. This occurs at many conventional fixed points of NLcrMs 
which are perturbatively accessible within the usual e-expansion in d = 2 + e dimensions. Since such 
operators are not prohibited from occurring in the action, they appear to threaten the very exis- 
tence of such fixed points. At the same time, for NLcrMs describing metal-insulator transitions of 
Anderson localization in electronic conductors, the strong RG-relevance of these operators has been 
previously related to statistical properties of the conductance of samples of large finite size ( "con- 
ductance fluctuations"). In this paper, we analyze this question, not for perturbative RG treatments 
of NLcrMs, but for two-dimensional Wess-Zumino-Witten (WZW) models at level k, perturbatively 
in the current-current interaction of the Noether current ( "non-Abelian Thirring/Gross-Neveu mod- 
els"). WZW models are special ("Principal Crural") NLcrMs on a Lie Group G with a WZW term at 
level k. In these models the role of high-gradient operators is played by homogeneous polynomials 
of order 2s in the Noether currents, whose scaling dimensions we analyze. For the Lie Supergroup 
G — GL(2iV|2iV) and k — 1, this corresponds to time-reversal invariant problems of Anderson local- 
ization in the so-called chiral symmetry classes, and the strength of the current-current interaction, 
a measure of the strength of disorder, is known to be completely marginal (for any k). We find that 
all high-gradient (polynomial) operators are, to one loop order, irrelevant or relevant depending on 
the sign of that interaction. 



I. INTRODUCTION 



Fluctuations are known to play a key role in suf- 
ficiently low-dimensional systems, whether classical or 
quantum, as they can preempt spontaneous symmetry 
breaking. When the symmetry is both global and contin- 
uous, the tool of choice to address the role of fluctuations 
in low-dimensional systems is the non-linear sigma model 
(NLcrM). However, the usefulness of NLcrMs has come 
to transcend situations in which a pattern of symmetry 
breaking is immediately obvious. For example, NLcrMs 
have been used with success in the context of Ander- 
son localization (see Ref. |l| for a review) to access the 
transition from a metallic to an insulating phase induced 
by weak disorder or to comp.utc probability distributions 
of spectral,0 wavefunction,0 and transport characteris- 
tics in chaotic metallic grains and disordered electronic 
systems B 

Quite generally, the construction of a generic NLcrM on 
a connected Riemannian manifold 371 of finite dimension 
n, the "target manifold", can proceed in the following 
way.u One assigns to any point from Euclidean space in 
d dimensions, specified by coordinates x^ (fj, = 1, ■ ■ ■ ,d), 
a point in the manifold 971 with the coordinates <fi l (x) 
(i = 1, • • • ,n). The simplest action S, which is made 
of two derivatives of the coordinates <ff , and is invari- 
ant under both the rotations of Euclidean space and 



reparametrization of the target manifold, is 
1 f d d a 



S = 



4nt 



..d-2 « 



G ij [ ( j ) (x)]d ll f(x)d^(x) (1.1) 



j [4>] is a component of the metric tensor on 971, 



where G 

t is the coupling constant, and a is the short-distance 
cutoff. 

The target manifold can be either compact or non- 
compact. An example of a NLcrM on a compact target 
manifold is the 0{N)/0(N - 1) NLcrM with 2 < N = 
3, 4, 5, • ■ ■ when the target manifold is the unit sphere 
gN-i j n TV-dimensional Euclidean space. When N = 3 
it describes spontaneous symmetry breaking in a classi- 
cal fcrromagnct. Non-compact target manifolds are of 
relevance to the problem of Anderson localization in the 
bosonic "replica limit" N — > (Lor when the manifold is 
generalized to a supermanifoldu In Anderson localization 
the coupling constant t has the meaning of the inverse of 
the mean dimcnsionlcss conductance^ 



The implicit assumption made in the construction (1.1) 
is that all the invariant scalars that contain 2s (1 < s = 
2, 3, ■ ■ • ) derivatives of the field can be ignored. The stan- 
dard justification for this assumption is that their "en- 
gineering dimension" 2s is much larger than the spatial 
dimension d = 2 + e, i.e., they are irrelevant in the renor- 
malization group (RG) sense, and this is expected to re- 
main so after renormalization in d = 2 + e dimensions for 
small e, and thus small t. 
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This assumption was called into question in Refs. 
[l6| , for which the main results can be illustrated most 
simply by the example of the 0(N)/0(N - 1) NLctM. 
We recall that the 0(N)/0(N - 1) NLerM has an infra- 
red unstable fixed point located, to one loop order, at 
t* = e/(N — 2), from which emerges a renormalization 
group (RG) flow to strong and weak coupling. In Rcf. [|, 
a family of perturbations of the 0(N)/0(N - 1) NLctM 
action (JO), which we shall call high-gradient operators, 
was considered. A high-gradient operator of order s is 
a homogeneous polynomial of order 2s in the derivatives 
of the fields (all located at the same point) which is a 
scalar with respect to both the symmetry group of the 
NLerM [i.e., O(iV)] and the rotation group of Euclidean 
space. The minimum (i.e., dominant, or "leading") value 
of the one-loop scaling dimcnsionsliZI of the high-gradient 
operators of order s at the fixed point t* is founda to be 



2s-s(s-l)t* +0{e 2 ). 



(1.2) 



Although strongly irrelevant by power counting (i.e., in 
the absence of fluctuation corrections, t* — > 0), high- 
gradient operators of order s thus acquire a one-loop 
scaling dimension smaller than two when the order 2s of 
derivatives is large enough so that st* ps sej (N — 2) ~ 2, 
and thus would appear to become relevant, based on the 
one-loop result. In d .=■ 2 dimensions, the lowest one- 
loop scaling dimensionta for all high-gradient operators 
of order s is 



2s-s(s- l)t + 0{t 2 ) 



(1.3) 



along the trajectory to strong coupling away from the 
infra-red unstable fixed point t = 0. Two-loop counter, 
parts to Eqs. fll.2| ) and (Oj) yield the same conclusion:Ej 
high-gradient operators of sufficiently high-order s ap- 
pear to be relevant for any given dimension d = 2 + e at 
the non-trivial fixed point. 

Similar results hold for the NLcrMs defined on the 
compact target manifolds (M and-JV are positive inte- 
gers) $tdM+ N)/Sp{M) x Sp(N)mJJ{M+N)/\J{M) x 
U(iV)|Jta 0(M + N)/0{M) x or/v ^^ pm H on families of 
compact Kahlcr (and super) manifolds £j Generalizations 
to the non-compact target manifolds Sp(M, N)/Sp(M) x 
Sp(JV), U(M,N)/XJ(M) x U(A), and 0(M,N)/0(M) x 
O(N) follow from the rule that the coupling t of the 
compact NLerM entering in one-loop anomalous dimen- 
sions must be replaced by — t in the corresponding non- 
compact NLctM. In Anderson localization, compact tar- 
get manifolds arise when using fermionic replicas for dis- 
order averaging, whereas non-compact target manifolds 
arise when using the bosonic replicas for disorder aver- 
aging. If one uses supersymmetric disorder averaging, 
the resulting NLerM has both compact and non-compact 
sectorsJd The high-gradient operators in the NLerM de- 
fined on AdSs x S 5 (AdSs is non-compact whereas S 5 is 
compact) have also been discussed in the context of the 
AdS/CFT correspondence. (See, for example, Refs. [T^ 

0) 



The substitution t — >• — t does not affect the value of 
the minimal (i.e., dominant, or "leading") one-loop scal- 
ing dimensions, _w hen the spectrum of anomalous one- 
loop dimensionsEJ of all high gradient operators of order 
s is distributed symmetrically about zero. This turns 
out to be the case whenever m, n > 1 in the above exam- 
ples. On the other hand, there are some target manifolds, 
the simplest examples being S 1 "" 1 = 0(N)/0(N - 1) 
and CP"" 1 = V(N)/V(N - 1) x U(l), for which the 
full spectrum of one-loop anomalous dimensions of order 
s turns out to be not symmetric about zero, in which 
case the substitution t — > — t matters. For example, 
high-gradient operators are made more irrelevant by one- 
loop renormalization effects in the non-compact NLerM 
on U(AT-1, l)/U(iV-l) x U(l). (We refer the reader to 
Appendix |A| for a more detailed discussion of "one-sided" 
versus "two-sided" spectra of one-loop anomalous scaling 
dimensions for high-gradient operators in NLcrMs.) 

Of course, one can only conclude that high-gradient 
operators become relevant for sufficiently large values of 
s, if the strong relevance seen in the one- loop expressions 
for their scaling dimensions persists when all higher loop 
contributions (not computed here or in other works on 
this subject) have been taking into account. For example, 
the one-loop expressions may not be characteristic in the 
large- s4i mi t, if the actual expansion parameter is not e 
; e £3'E2l As any insight for resolving the nature of the 



but sei 

e expansion for high-gradient operators in NLcrMs must 
come from outside the e expansion itself, progress has 
stalled since the early 1990's. 

The aim of this paper is to study the operators that 
play the role of the high-gradient operators in field the- 
ories which are|—twxLjdimensional Wess-Zumino-Wittcn 
(WZW) theoriescj l3 on a Lie group G, perturbed by an 
interaction quadratic in the Noether currents ("current- 
current interaction"). Such theories are often also re- 
ferred to as "two-dimensional non-Abelian Thirring (or 
Gross-Neveu) models". Any WZW theory, which is a 
Principal-Chiral-Non-Lincar-sigma model supplemented 
by a WZW term at level k, gives a prescription to 
construct high-gradient operators in terms of powers of 
Noether currents. Because it is possible to represent 
the Noether currents in WZW theories in terms of free 
fermions,c3 one might be inclined to think that such op- 
erators are perhaps not capable of displaying a "patho- 
logical" spectrum of scaling dimension as in Eq. ( p_.2j ). 
However, as we demonstrate in this paper, the situation 
is more interesting. Indeed, we will see that under con- 
ditions specified below, the one-loop spectra of the form 
(1.3) and ( [L~2] ) can be realized by perturbing a WZW 
critical point by a current-current perturbation. 

We also want to investigate if there is a difference be- 
tween the properties of high-gradient operators in unitary 
and non-unitary non-Abelian Thirring models. This is 
important because NLcrMs describing the physics of An- 
derson localization are non-unitary field theories. More- 
over, high-gradient operators in these theories have been 
previously related to the statistical fluctuations of the 
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conductance of a disordered metal.@ _ ii In this con- 
text, an appealing physical interpretation of the spectra 
(1.3) and ( [T^ ) has been proposed, attributing them to 
a broad tail in the probability distribution of the con- 
ductance. Ej EHI However, given that this interpretation 
depends crucially on the ability to invert the s — > oo and 
e — > lim its in spectra which are analogous to those in 
Eqs. ( |l.3| ) an d (P--2|), it would be useful to have an ex- 
ample of a critical field theory describing a problem of 
Anderson localization for which one can study the RG- 
relevance of high-gradient operators without resorting to 
the e-expansion, and for which one can reasonably expect 
a broad distribution of the conductance. 

We now provide an outline of the article and a sum- 
mary of our results. 

It is shown in Sec. |l| that high-gradient operators in 
the (unitary) su(2) fc Thirring model with strength g of 
the "current-current interaction" , are made more irrele- 
vant by the presence of these interactions when the latter 
are (marginally) irrelevant (g < 0, in our conventions). 
On the other hand, along the renormalization group (RG) 
flow driven by a (marginally) relevant current-current in- 
teraction (g > 0, in our conventions), a one-loop spec- 
trum of the form (1.3) is recovered in the "classical" limit 
1/k — > 0. The inverse level 1/k plays here the role of a 
"quantum" parameter. Indeed, for any finite k, we find 
that the quadratic growth in s in the unbounded one- 
loop spectrum (1.3) does not persist for values of s larger 



than k. In effect, 1/k determines the efficiency in "tam- 
ing" the strong RG-relevance of high-gradient operators 
seen at one-loop order, which is related to the fact that 
there exists a representation of the current algebra of the 
level- k WZW theory in terms of free fermions. 

Section [II is devoted to high-gradient operators in 
what we will call the g\(M\M) k Thirring (or Gross- 
Neveu) model which was discussed in Rcf. |33j. This 
is the g\{M\M) k WZW theory on the Lie Supergroup 
GL(A/|Af), perturbed by two current-current perturba- 
tions, one which we call <? M , which is exactly marginal, 
and another which we call g A , which flows logarithmi- 
cally under the RG at a rate dependent on g M . In spite 
of the presence of an RG flow of the coupling g A there 
exists a sector of the theory, the so-called PSL(M|M) 
sector, which is scale (conformally) invariant through- 
out.!^ The high-gradient operators turn out to reside in 
this conformally invariant sector, and are unaffected by 
the presence of the coupling g A . We will show that, for 
k = 1, the spectrum of one-loop anomalous dimensions 
of high-gradient operators is fundamentally different for 
positive and negative values of the coupling constant g M . 
In particular, when g M > all high-gradient operators 
arc made more irrelevant by the current-current per- 
turbations, whereas they are made more relevant when 
g M < 0. We close Sec. [II by comparing the anoma- 
lous scaling dimensions of high-gradient operators in the 
g\(M\M) k Thirring (or Gross-Nevcu) models and those 
in the GL(2A|2A)/OSp(2A|2A) NLcrMs, observing that 
they behave in the same way. 



The result that the spectrum of one-loop anomalous 
dimensions of high-gradient operators is strongly depen- 
dent on the sign of g M has important implications in the 
context of Anderson localization because, as discussed in 
Ref. H, the g\{M\M) k Thirring (or Gross-Nevcu) model 
at k = 1 describes a disordered electronic system, where 
g A > and g M > correspond to the strengths of dis- 
order potentials. The theory with g M > thus offers 
an example of a critical theory for Anderson localization 
with no relevant high-gradient operator. - For example, 
this field theory describesE3 a tight-binding model of elec- 
trons on the honeycomb lattice with (real-valued) ran- 
dom hopping matrix elements which are non-vanishing 
only between the two sublattices of the bipartite hon- 
eycomb lattice (see also Ref. [}4|). Versions of the hon- 
eycomb tight-binding model provide the basic electronic 
structure of graphene. In the classification scheme of 
Zirnbauer, and Altland and Zirnbaucr.Bj E3 this model 
belongs to the "chiral-orthogonal" symmetry class (class 
BDI). (Another example of a problem of Anderson lo- 
calization in the same symmetry class is provided by a 
random tight-binding model .ati a square lattice with ir- 
flux through every plaquetteo) 

By contrast, when g M < 0, the spectrum of one-loop 
scaling dimensions is unbo unded from below for any k 
as is the case in Eq. (1.2). The full spectrum of one- 
loop anomalous scaling dimensions of high gradient oper- 
ators as it appears, e.g., in the Grassmanian NLcrMs with 
target manifolds Sp(M + A)/Sp(M) x Sp(iV), U(M + 
N)/U(M) x U(A), 0(M + N)/0(M) x O(A), which is 
symmetric about zero, is only recovered in the extreme 
"classical" limit M, k — > oo. In the context of Anderson 
localization, the case with g M < describes the surface 
state of a three-dimensional topological insulator in the 
chiral-symplcctic class (symmetry class CII) of Anderson 
localizationM9c3 

After concluding in Sec. |Tv| , we review in Appendix [B] 
the realization of the gl(27V|2A) fc=1 Thirring (or Gross- 
Neveu) model as a problem of Anderson localization in 
two dimensions in symmetry class BDI, which was estab- 
lished in Ref. Isl 



II. HIGH-GRADIENT OPERATORS AND su(2) fc 
WZW THEORIES 



The 0(3)/0(2) NLcrM with coupling constant t is the 
simplest example of a NLcrM containing infinitely many 
high-gradient operators all of which would appear to be- 
come relevant based on one-loop results. This happens at 
the infra-red unstable fixed point t*=eind = 2 + e>2 
dimensions within the one-loop approximation as long 
as the order s of these high-gradient operators is large 
enough. A precursor to this perturbative property also 
occurs in d = 2 dimensions close to the infra-red unstable 
fixed point t = as the NLcrM flows to strong coupling- 
Along this flow, the spectrum of one-loop dimensionstSI 
for the high-gradient operators is unbounded from below. 
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In two dimensions, the 0(3)/0(2) NLcrM, supple- 
mented by a topological theta-term at 6 = 7r, flows 
to a critical field theory, the SU(2) WZW theory.jp 
su(2) fc=1 current algebra, su(2) fc=1 WZW theory.a - 
The strongly relevant high-gradient operators near the 
infra-red unstable fixed point t = must become irrel- 
evant at the WZW critical point, because the full op- 
erator content of the su(2) fe=1 WZW theory is known 
to contain only a finite number of relevant fields (with 
scaling dimensions bounded from below and above by 
zero and two, respectively). The purpose of this section 
is to perturb the su(2) fe WZW theory with a current- 
current perturbation and to examine the fate of those 
operators in the su(2) fc WZW theory which correspond 
to the high-gradient operators in the 0(3)/0(2) NLcrM. 
We will refer to these operators still as "high-gradient 
operators" . We are going to argue that the spectrum 
of one-loop scaling dimensions associated with all high- 
gradient operators is bounded from below by the low- 
est one-loop scaling dimension corresponding to high- 
gradient operators of order k. This result is very dif- 
ferent from the unbounded spectrum (1.3) of one- loop 



scaling dimensions associated with high-gradient opera- 
tors in the two-dimensional 0(3)/0(2) NLcrM. 

In the following sections, we first review the su(2) fc 
WZW theory perturbed by a current-current interaction. 
Second, we identify high-gradient operators of order s. 
Finally, we compute the leading one-loop dimensions of 
high-gradient operators of order s up to one loop. 



A. Definitions 

The most fundamental property of the su(2) fe WZW 
theory is the existence of a pair of holomorphic and anti- 
holomorphic SU(2) Noether currents, J x , J 2 , J 3 , and J 1; 
J 2 , J$, respectively, which satisfy the affine (Kac- Moody) 
current algebra 



kC. 



J a (z)Jp(0)=0, 



kC„o i 



1 + y af ruo)+ 



(2.1a) 



at level k = 1,2,3, ••• , where the invariant (Casimir) 
tensor of rank 2 in su(2) has the contravariant and co- 
variant representations (in our conventions) 



C. 



a/3 



1, 

2 aP ' 



C afS = 25 



(2.1b) 



respectively, while the structure constant of su(2) is the 
fully antisymmetric Levi-Civita tensor of rank 3, 



f 7 
Jap 



£ q/3 7 i a : &7 = 1,2,3. 



(2.1c) 



The dots in Eq. ( |2.1a ) are terms of order zero and higher 
in powers of z (z) with z = x + iy (z — x — iy) the holo- 
morphic (antiholomorphic) coordinates of the Euclidean 



plane. We shall also refer to the (anti) holomorphic sec- 
tor of the theory as the (right- ) lef t-moving sector. 

The su(2) fc current algebra ( |2.l[ ) has a representation 
in terms of free-fermions. More precisely, it is obtained 
from the actiorc3 



s * : =E 



^(^Mci+tldtJ (2-2a) 



constructed from fc-independent flavors of left (?/>) 
and right ($) moving Dirac fcrmions, whereby each 
one transforms in the fundamental representation of 
SU(2) xSU(fc), with the partition function 



Z, 



V[nl>\l>,i>\$\ exp(-SJ. 



(2.2b) 



One has the operator product expansions (OPE) 



i(^)^(o) = f;Wi(o) 



S.,,6. 



i' cd 



i(# d J(o)=^W^(o) 



z 

S.A 



U cd 



(2.3) 



<U*)#J(o)~o, 



for 1, d = 1, • • • , k and c, d = 1, 2. In turn, the OPE ( gj ) 
imply that the left and right Noether currents 



t=l " 1=1 

(2.4) 



with a = 1,2,3 obey the SU(2) fe current algebra (2.1). 
The field theory defined by Eq. (2.2) is a frcc-fcrmion 



field theory. The content of local operators is thus known. 
It contains a finite number of fields whose scaling dimen- 
sions are bounded between and 2 and are thus relevant, 
as it should be for a field theory defined on a Hilbert space 
with a positive definite inner product and with a spec- 
trum bounded from below which is built on the Dirac- 
Fermi sea, in short a unitary field theory. Clearly, within 
the set of powers of the Noether currents (|]J) there is 
thus no room for an infinite family of relevant operators. 

We perturb the free-fermion field theory by a current- 
current interaction O i of the SU(2) fe currents (2.4). 



Z:= / V[4>\$,$\$\ exp(-S). 
dzdz 



S '■= + g 



2iri 



■Oi(z,z), 



(2.5) 



Oj{z,z) := C a »J a (z)J p (z) = 2J a {z)J a {z). 



We take the coupling constant g to be real. The (uni- 
tary) field theory (^3) is often referred to as a 11011- 
Abelian Thirring ( Gross- Neveu) model. Suit able non- 
unitary generalizations of the field theory (2.5) compute 
(disorder average) moments of Green's functions in a 
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class of problems of Anderson localization in d = 2 di- 
mensions that we will investigate later on in this paper. 
The one-loop beta function, 



8 =^ = 
1 g dl 



(2.6) 



encodes the change in the coupling constant caused by 
the infinitesimal rescaling a — > (1 + dl)a of the short- 
distance cutoff a. Thus, the current-current interaction 
is (marginally) irrelevant (relevant) for g < (g > 0) 
with the free-fermion fixed point at g = 0. 

The SU(2) fc Noether currents (2.4) are those appear- 
ing at the non-trivial fixed point of the Principal Chiral 
NLaM-oa-the SU(2) group manifold with a Wess-Zumino 
term£3 _ EZI This has, the well-known (Euclidean) action 



k 
16V 



J d 2 xtr (d^G^d^G) + kT[G], (2.7a) 



where G € SU(2) is a group element, and the integral 
T[G] over a three-dimensional ball B with coordinates r 



and whose boundary dB is d ■■ 
space, 



2-dimensional Euclidean 



m : =^J ^re^ti (G~ 1 d fl GG~~ 1 d l ,GG~ 1 d x G) 

□ (2.7b) 
is the Wess-Zumino term.Ej 

The Noether currents which generate the 
SU(2) left x SU(2) right symmetry at the critical point 
of the WZW theory can be ful ly represented by the 
fermionic expressions in Eq. (2.4). In the bosonic (i.e., 
NLerM) representation, these currents are built out of 
first-order derivatives of the bosonic fields, 

J a oc fctr [{dG)G- x a a ] , J a oc fctr [G- 1 (dG)a a ] . 

_ (2-8) 
The relationship (Bj|) suggests that comp osite opera- 
tors built out of monomials in the currents (2.4) in the 
WZW theory are the counterparts of the high-gradient 
operators in NLctM. For this reason, we shall still call the 
former family of composite operators high-gradient oper- 
ators. The "classical" counterparts of the high-gradient 
operators of order s in the NLctM are thus the homoge- 
neous polynomials 



rpo 



■Qft, ■■■Ct s 



J ™ ' ' ' J r\. 



J, 



(2-9) 



of the left and right currents that arc invariant under the 
diagonal SU(2) symmetry group of the interacting the- 
ory.c3 The generating set of classical high-gradient oper- 
ators of order s is specified once all the linearly indepen- 
dent rank 2s tensors T a ^"" yS '" in the adjoint representa- 
tion of SU(2) that are invariant under SU(2) transforma- 
tions can be fully enumerated. In turn, the most general 
SU(2) invariant tensor of even rank in the adjoint rep- 
resentation is the product of the Casimir tensor of rank 
2H 



The high-gradient operators in Eq. ( |2.9| ) are classical in 
the sense that quantum fluctuations encoded through the 
Pauli principle (or, equivalently, through the underl ying 
Dirac-Fermi sea) in the free-fermion representation (2.4) 
of the current algebra, have not yet been accounted for. 
To account for these quantum fluctuations, one needs to 
introduce a point-splitting procedure that allows for the 
proper normal ordering, i.e., JJae correct subtraction of 
all short-distance singularitiesuj 



lim 



rpC 



(all short-distance singularities) 



(2-10) 



Two objects of the form (2.9) that are linearly inde- 
pendent classically might not survive as a pair of dis- 
tinct quantum operators of the form ( [2.1C ) after nor- 
mal ordering has been implemented. More precisely, one 
might anticipate that the underlying free-fermion repre- 
sentation of the current algebra must manifest itself as 
soon as the order s becomes larger than the number k of 
fermionic flavors by changing the book-keeping relating 
classical expressions labeled by SU(2) tensors of rank 2s 
and quantum operators. 

Indeed, we are going to show that this is the mecha- 
nism that prevents high-gradient operators of order s > k 
from acquiring one-loop scaling dimensions smaller than 
the smallest one-loop scaling dimensions associated with 
the set of all high-gradient operators of order s < k. 
In other words, the smallest one-loop dimension associ- 
ated with the set of all high-gradient operators is reached 
within the set of all high-gradient operators of order 
s < k when g > 0. It is thus bounded from below when 
g>0. 

Had we ignored the underlying free-fermion represen- 
tation of the current algebra altogether, we would have 
wrongly predicted that, when g > 0, the one- loop di- 
mensions associated with the classica l ob jects (2.9) are 
of a form similar to the ones in Eq. (1.2) i.e., that the 
set of one-loop dimensions of high-gradient operators is 
unbounded from below. On the other hand, this classical 
prediction is recovered in the limit k — > oo with s/k — > 0. 
For this reason we shall separate the computation of the 
most relevant one-loop dimension associated with high- 
gradient operators of order s into the case when s < k 
and the case when k < s. 

In this context, we would like to remind the reader that 
the su(2) fc WZW theories are known to describe quantum 
critical points in the parameter space of quantum spin-S* 
antiferromagnetic chains when k = 2S. Here, we observe 
that both, the number of relevant perturbations and the 
number of independent local composite operators built 
out of the generators of SU(2) which are SU(2) singlets, 
grows with S. [For 5 = 1/2 the algebra obeyed by the 
Pauli matrices only allows one invariant SU(2) tensor of 
rank 2, the 2x2 unit matrix.] On the other hand, the 
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strength of quantum fluctuations in su(2) fc WZW theo- 
ries decreases with increasing k = 2S for the same rea- 
son as the role of quantum fluctuations decreases with 
increasing S in quantum spin chains. 



B. Anomalous dimensions of high-gradient 
operators 

As the most general SU(2) invariant tensor of even 
rank in the adjoint representation is the product of the 
Casimir tensor of rank 2,e3 we define the three diagonal 
SU(2) invariants out of the three current bilinears 

H := C a/3 J a Jp, A := C a ^J a Jp, B := C Q/3 J a Jp, 

(2.11a) 



together with the SU(2) invariant 
C := AB. 



(2.11b) 



The space of the high-gradient operators is then spanned 
by the family 

Sh\H'- 2 C,--- ,H 2 C [s/2] -\C [s/2] } (2.12) 



made of [s/2] + l "classical" operators.^] We call these op- 
erators "classical" because we have not yet taken into ac- 
count the short distance singularities associated with the 
definition of composite operators (i.e., the "Pauli princi- 
ple" discussed above). As announced below Eq. ( [2.10 ), 
these singularities need to be subtracted from the "clas- 
sical" expressions (2.12) upon normal ordering. We shall 



nevertheless ignore the issue of normal ordering at first 
and compute the one-loop RG equation for these "un- 
rcgularized" (or un-normal-ordered) operators, a step of 
no consequence in the ("classical") limit s/k — > 0. We 
shall then contrast this un-regularized calculation with 
the full quantum calculation for the special case oik = 1 , 
i.e., when the proper normal ordering procedure has been 
accounted for. 

We shall see that the calculation without normal or- 
dering gives an infinite tower of high-gradient operators 
that are all relevant to one-loop order for sufficiently large 
s and for g > 0. The one- loop spectrum of anomalous 
dimensions is identical to the one in the 0(N)/0(N — 1) 
NLerM when N = 3. Indeed, once the normal or dering 
procedure is ignored, the high-gradient operators (2.12) 
are analogous to those discussed in Ref. [15[ and the calcu- 
lations of the anomalous dimensions in the su(2) fe WZW 
model and in the 0(N)/0(N - 1) NLcrM run along par- 
allel tracks. 

The effect of normal ordering is weaker the smaller s/k 
is, i.e., the closer proximity to the semi-classical limit of 
the WZW theory. To see this, consider the case when 
k > s. The "classical" expression T a i'" a * J ■ ■ - J s (z, z) 
for the composite operator made of a local product of 
holomorphic and antiholomorphic currents is modified 



upon normal ordering. To leading order in a short- 
distance expansion, this classical expression is replaced 
by 

k 

E 7 " , t .:(*) 

x ■■■)„ , ••••/„ ,. :(*) • ••• • (2-13) 

Here, the terms included in the • • • arise from the OPE 
for the product J {z)J (0) when any two flavor in- 

i i j j 

dices i i and i • are identical. Evidently, normal ordering 
(or the Pauli principle) has a much more potent effect 
when k < s, for the condition l 1 ^ • • • ^ i s ^ i x ^ ■ • • ^ 
l s can then never be met so that the leading order term 
above is absent. The operator contents with and without 
normal ordering thus look very different. When k < s, 
some operators in the set ( 2.12| ) completely disappear to 
leading order because of Fermi statistics. This will be 
demo nstrate d ex plicitly for the case of k = 1 [see Eqs. 
(2.28) and (2.29) below], for which we will show, after 
correctly taking into account normal ordering, that all 
high-gradient operators which would be relevant classi- 
cally (when g > 0) disappear from the operator content. 



1. RG equation for un-regularized high-gradient operators 



To compute the leading one- loop scaling dimensions 
for the high- grad ient operators ( ^.12| ), we start from the 
field theory ( |2.5| ) in which we substitute the action by 



S '. — S L 



dzdz 
2m 



2m-\-n—s 



2s-2 



m,n— 



dzdz 
2ni 



(2.14) 



To determine the one-loop dimensions of the couplings 
{Zm,n|2m + n = s}, we do not need the full one-loop RG 
flows, i.e., the RG equations for the coupling constants 

(s) (s) 

up to and including order Zm'nZp,q, but only the linear 



(s) 

in Zffi ji contributions to the one-loop RG flows. Thus, all 
we need arc the OPE of C m H n {z, z) with 0,(0), where 
the integers m and n satisfy 1 < 2m + n = s < k. Fur- 
thermore, we shall introduce the short-hand notation 

AxB = C^> A(z, z)B(0) = —_C(z, !) + ••■ (2.15) 

zz 

for the OPE relating the operators A, B, and C. Here, 
the dots are meant to contain not only regular terms of 
zcroth and higher order in z or z but also second and 
higher order poles in z or z. 

As an intermediary step, one verifies that the OPE 
( p. 15 ) between the building blocks H and C with Oj 
(observe that Oj = H) are 



HxOj 



-AH, CxOj = 0. 



(2.16) 
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Here, we introduced yet another short-hand notation 
i 1 

A ■ ■ ■ B or A ■ • • B, by which we mean that one current 
in A and one current in B arc contracted with the rule 



C, 



c 



(2.17) 



for any a,/3 = 1,2,3 and i, t! = 1, ■•• ,k at the free- 
fermion fixed point g = 0, where 



= r i(^al^ J Qt:= ^t(^k!^. (2.18) 



When A and B consist of more than one J a or J a , and 

when there arc many possible Wick contractions between 

■ 1 

A and £>, the short-hand notations A - ■ ■ B, A - ■ ■ B and 



■ 1 



A ■ ■ ■ B mean the resulting operator obtained by taking all 

possible s uch W ick contractions. One also verifies that 
the OPE fl2.15D between the building blocks HH, CH, 
and CC with Oj are 



HHxOj = -AH 2 +4C*, 



—I t~ 



(2.19) 



CHxOj = CCxOj=0. 



We then infer that, for any pair (to, n) of positive integer 
that satisfies 1 < 2m + n = s < k, 



C m H n xO I = CxC?! xmC m - 1 H n 



HxOjX nC m H n - 1 



rn — 1 Tjn — 1 



CHxO I x mnC m - L H 



+ GCxOj x 



m(m-l) cm _ 2Hn (2.20) 



+ HHxOj x n(n V iT- 2 

= -2n {n + l)C m H n 
+ 2n{n-l)C m+1 H n - 2 . 

The contributions to the RG equations obeyed by the 
couplings Zra, n where 1 < 2m + n = s < k needed to 
extract the spectrum of one-loop dimensions are 

(2-2*)zW n + 4gn(n + l)zW n 



dl 



(2.21) 



4 ff (n + 2)(n + l)Z« +a + 



Here, the dots include non-linear contributions of second 
order in g or Zm] n . 

The linearized RG flows (2.21) are closed. This is a 
justification a posteriori for neglecting the RG effects of 



current monomials with repeating flavor indices. The 
linearized RG flows ( 2.21 ) have a lower triangular struc- 
ture, i.e., there is no feedback effect on the flow of a 
high-gradient operator of order s from lower-order high- 
gradient operators. Thus, we conclude that the leading 
[s/2] + 1 one- loop scaling dimensio ns ass ociated with the 
family of high-gradient operators ( 2.12} ) when k > s = 
2m + n are given by 



2(2m + n) - 4gn(n + 1). (2.22) 



Observe that the spectrum of anomalous dimensio 



/m,n 



-Agn(n + 1), 



2m 



ns0 

(2.23) 



is one sided with respect to 0. When g < these anoma- 
lous dimensions are positive, i.e., the scaling dimensions 
are larger than their engineering value. The opposite 
happens when g > 0, i.e., when the current-current per- 
turbation is (marginally) relevant. When g > and for 
a given 1 < s < k, the smallest one-loop anomalous di- 
mension occurs for the pair (m, n) = (0, s), 



W 

^min 



mm 

2m+n— s 



?s(s + l). (2.24) 



For g > 0, the quadratic dependence on s can overcome 
the linear dependence on s in the one-loop dimension 
x min : ~ 2s + 7^j n . If the order s (1 < s < k) is allowed 

(s) 

to be sufficiently large, the one-loop dimension x y m ( n de- 
creases past the value 2 and eventually becomes negative. 
The quadratic dependence on s is reminiscent of that 
for the one-loop dimensions (1.2) in the 0(N)/0(N — 1) 



NLerM. However, in contrast to the (2 + e)-dimensional 
0(N)/0(N-1) NLerM at its non-trivial fixed point t*, a 
value smaller than 2 for the one-loop dimensions a^; n 
is not a threat to the internal stability of the WZW 
fixed point g = since it occurs along a flow to strong 
coupling. Moreover, it is known that in d = 2 dimen- 
sions the 0(3)/0(2) NLerM with thcta term at 6 = ir 
flows in the infrared into the level k = 1 SU(2) WZW 
fixed point. While the spectrum of one- loop dimensions 
of high-gradient operators at the WZW fixed point is 
bounded from below (as we will recall below), the spec- 
trum of these operators is unbounded from below in 
the weakly coupled 2-dimensional 0(3)/0(2) NLerM (the 
presence of the theta term docs not affect this result). 



2. Normal ordering revisited 

We shall illustrate the effects of the Fermi statistics for 
the family of high-gradient operators ( 2.12 ) when s = 2 
for the case of a (marginally) relevant (g > 0) current- 
current interaction. We shall then show for the special 
case of k = 1 and s = 2 that the two one-loop dimen- 
sions associated with the family of high-gradient opera- 
tors (2.12) arc unchanged, to one loop order, i.e., 



C 0,2 



(*) 
c l,0 



(2.25) 
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We sta rt from the family of high-gradient opera- 
tors ( 2.12 ) with s = 2. For clarity of presentation, we 
rename the two members of this family, 

1 = C a ^ J a JpC lS J^Jfi, 2 = C a ^ J a J pC lS J^Jfr. 

(2.26) 

As implied by Eq. ( 2.10 ) these are two classical expres- 
sions. The two quantum expression s invo lve point split- 
ting and normal ordering as in Eq. ( [2. 10 ). 

Without loss of generality, we consider only the left 
current sector. Normal ordering of 

kC„ 



Zq 2 ; associated with : 2 :, we find that their one- loop 
dimensions remain equal to their engineering dimensions, 



J Q (z)J^(0) + -e^ 7 J 7 (0) + -e Q07 aJ 7 (O) 



i.=i 



(2) 
"1,0 



J 2 ) 



(2.31) 



The lesson that we draw from the example s = 2 
and k = 1 is that it is necessary to use normal order- 
ing to properly define composite operators. Had we not 
used normal ordering, we would have incorrectly pre- 
dicted that there are infinitely many high-gradient oper- 
ators which become relevant, at one-loop order, for large 
enough s and for g > 0. We believe that for a generic 
value of k, there is no infinity of one- loop relevant high- 
gradient operators. Only a finite number of high-gradient 
operators become relevant, at one-loop order, for large 
enough s and for g > when k > 1. 

In the next section, we turn attention to a non-unitary 
WZW model of relevance to the problem of Anderson 
localization to investigate whether the loss of unitarity 
opens the door to an infinity of relevant high-gradient 
operators. 



(2.27) 



amounts to the subtraction from Eq. ( 2.27 ) of the terms 
singular in the limit z — > 0, 



+ -f E : {vrHa. - - rH a ^ b M ■ (o) 
i=i 

(2.28) 



for a, ft = 1, 2, 3. The prop er quantum interpretation of 
the classical currents ( 2.26 ) is then 



O i: (M) = 4 E •J a J a :{z):JpJ fi :(z), 

<*,/3=l 

3 (2.29) 
2 :(z,z) = 4 E : J a J : (z) : Jj p : (z). 

a, = 1 



3. High-gradient operators when k = 1 

When k = 1, the summation over unequal flavors dis- 
appears in Eq. (2.28). (Observe in passing that : 1 : 
is then proportional to one component of the energy- 
momentum stress tensor.) One then verifies the OPE 



:0 1 : xO I =3:0 1 :-9:0 2 :, 
:0 2 : xO I =:0 1 :-3:0 2 : . 



(2.30) 



If we diagonalize the linearized one-loop RG flows for 

(2) 

the coupling Z{ associated with : 1 : and the coupling 



III. HIGH-GRADIENT OPERATORS AND 

g\{M\M) k WZW THEORIES 

An interesting example of a problem of Anderson 
localization in two dimensions which possesses a spe- 
cial so-called sublattice (or chiral) symmetry (SLS) and 
TRS (thus belonging to the "chiral-orthogonal" symme- 
try class BDI in the classification scheme of Zirnbauer, 
and Altland and Zirnbauem^Lj) is as follows. Consider 
a tight-binding model of fcrmions on a honeycomb lat- 
tice with random real- valued hopping matrix elements 
of non- vanishing mean, which do not connect the same 
sublattice (so that SLS is preserved). [A related real- 
ization of the same problem of Anderson localization is 
provided by a random tight-binding model on a square 
lattice with flux-7r through every plaquette, introduced 
in Ref. [39[] In the absence of disorder this band struc- 
ture is known to exhibit the energy-momentum disper- 
sion law of two species of (relativistic) Dirac fermions at 
two points in the Brillouin zone at low energy near the 
Fermi level (at zero energy) . It was shown in Ref. |3^ that 
the SLS-prcscrving disorder discussed above leads to a 
theory for the disorder averages which, in the supersym- 
metric formulationjj is a GL(2.ZV|2A0 Thirring (Gross- 
Neveu) model. In other words, the problem of two- 
dimensional Anderson localization on the honeycomb lat- 
tice preserving SLS and TRS, is described by a set of 
Dirac fcrmions (and SUSY boson partners) perturbed by 
a current-current interaction of the Noether currents of 
its underlying GL(2A^|2A r ) (super) symmetry. The inter- 
action strength corresponds to the strength of the disor- 
der. The system of free P4p\c fermions (and SUSY boson 
partners) is well knownOc3 to be described by a WZW 
model on the supergroup Gh{2N\2N) with gl(2iV|2Ar) fc 
conformal Kac-Moody current algebra symmetry at level 
k = 1. 
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This section is devoted to the one-loop RG analysis 
of high-gradient operators in the perturbed gl(M\M)k 
WZW theory. The main result of this section and of this 
article applies to the case of level k = 1 of relevance for 
the random tight-binding models discussed above. In or- 
der to state this result, we first need to recall from Ref. 
||that the GL(2N\2N) Thirring (Gross-Neveu) models 
possess two coupling constants; one, g M , which does not 
flow under the renormalization group (RG) and another, 
g A , which flows logarithmically under the RG and a rate 
dependent on g M . Our main result then suggests that all 
higher-order gradient operators are more irrelevant in the 
presence of the current-current interaction with g M > 
than at zep, coupling g M = 0. A positive g M can be 
intcrpretcdEj as the variance of the disorder strength in 
the random tight-binding model in symmetry class BDI. 
For the opposite sign of the coupling constant g M < 0, 
on the other hand, higher-order gradient operators have 
a spectrum of one-loop dimension s th at is unbounded 
from below very much as in Eq. (1.2). In the context 
of Anderson localization, the case with g M < describes 
the surface state of a three-dimensional topological insu- 
lator in the chiral-symplectic. class (symmetry class CII) 
of Anderson localizationJlJ'liil 

As in Sec. [n| we are going to distinguish two limits. In 
the first (classical) limit, 



M oo, 



k 



oo, 



(3.1) 



OPEs between the high-gradient operators can be ob- 
tained without any reference to the composite nature of 
the currents. One then recovers a spectrum of one-loop 
scaling dimensions for high-gradient operators that mim- 
ics closely that of the NLerMs discussed above. The sec- 
ond limit, 



M = 1,2,3, 



1. 



(3.2) 



is the opposite extreme to the first one in that the nor- 
mal ordering of the currents and thus of the high-gradient 
operators is essential and changes dramatically the spec- 
trum of one-loop scaling dimensions from the "classical" 
limit (O). 



A. Definitions 



Our starting point is a two-dimensional coeformal field 
theory characterized by the current algcbraE3 



J A B (2)J C » 



BD 
AC 



+ L[^jn m + t BUjB m 



Jf (z) J C D (0) =^f- + \ [ifijf (0) + J* (0)] 



JaWc (0) =0, 



(3.3a) 



where 



and 



CIS := H 



9g=-H 



BCzB 
°Gi 



BD _ ^_^BC+D(1 



"AC 



(3.3b) 



(3.3c) 



with the indices A,B,C,D = 1, • • • ,M+N, where 5^ de- 
notes the Kronecker delta. The capitalized indices A, B, 
C, and D also carry a grade which is either for M out 
of the M + N values that they take or 1 for the remaining 
N values. It is the grade of the indices A and B that en- 
ters expressions such as (— ) A or (— ) AB . The grade (1) 
will shortly be associated with bosons (fermions). The 
positive integer k is t he le vel of the current algebra Q3.3| ) . 
The current algebra (3.3) is associated with the Lie su- 
pera lgebr a gl(M|A) defined by the structure constants 
Eq. (pic]) for A, B, C, D = !_, ■■■ ,M + N. When N = 0, 



the structure constants ( 3.3c ) reduce to 



9g 



*AC 



(3.4) 



for A,B,C,D = 1, • • • , M. These define the Lie algebra 
gl(M) of the non-compact Lie gr oup GL(M). When M = 
0, the structure constants ( |3.3cj ) reduce to 



Z AG 



°Ai 



(3.5) 



for A,B,C,D = 1, • ■ • ,N. These define the Lie algebra 
u(N) of the compact Lie group U(A). 

There exists a frcc-fermion and free-boson realisation 



of the current algebra (3.3) defined by the actionES 

dzdz 



5» 



E 

i=i 



2?ri 



{^ A \d^ AL + i> A \di> AL ) (3.6a) 



with the partition function 



Z„ := / V,V^]exp(-S J, 



(3.6b) 



where it is understood that ip Al and tp Al are complex- 
valued integration variables for the M values of A with 
grade while ip Al and ip AL are Grassmann- valued integra- 
tion variables for the N values of A with grade 1 , regard- 
less of the value take n by the flavor index i = 1, • • • , k. 
The current algebra (3.3) is then realized by the repre- 
sentation 



h 

E 

1=1 



i=i 

as follows from the OPE 

tMtfHo) = (-i) AB+ Vf (*)^(o) 



(3.8) 
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with A, B 



,M + N and i = 1, 



,k. 



The expressions in Eq. (3.7) form a representation of 
the g\(M\N) k current algebra in terms of free fermions. 
There are two Casimir invariants of rank 2 in g\(M\N) 
that we use to perturb the free field theory ( |3 .6[ ) with 
two types of current-current interactions, both of wJaich 
are invariant under the global GL(M\N) symmetry,E3 



Z 



S : — S, 

O a := 
O m := 



dzdz /g A 
2ni \~2tt 

A 7B 



O 



, 9m q 
2tt 



M 



(3.9) 



-^(-ir^(-i) 



-str Jstr J, 



■JaJb(-iY 



-str ( J J) 



Formally, one may allow the coupling constants g A and 
<7 M to take on any real (i.e., positive or negative) values. 
However, to make connection with the above mentioned 
two-dimensional tight-binding models in symmetry class 
BDI of Anderson localization, we must demand that g A 
and g M be positive. (See Appendix [b].) 

The "classical" counte rpar ts to the high-gradient op- 
erators of order s in Eq. (2.9) are the homogeneous poly- 
nomials 



■■A s Aj 

■B.B, 



-A, J B 1 
■B. J A, 



. . . J B ' T Bl 



(3.10) 



of the left and right currents that are invariant under the 
diagonal GL(M | N) symmetry group of the interacting 
theory.EZI The set of "classical" high-gradient operators 
of order s is specified once all the linearly independent 
rank 2s invariant tensors T^ 1 g 3 ^ 1 ^" in the adjoint 

representation of GL(M\N) which are invariant under 
GL(M\N) transformations have been enumerated. At 
the quantum level, normal ordering defines the quantum 
high-gradient operators of order s as in Eq. (2.1C). 

We are now going to specialize to the case M = N 
where the beta function for the coupling constant g M 
vanishes identically, an exact result Ej (As already men- 
tioned, the other coupling constant g A flows logarithmi- 
cally at a rate set by g M .) The sector which we loosely 
denote by 



PSL(M|M) - GL(M|M)/U(1) x U(l) 



(3.11) 



remains scale (conformally) invariant for any value of g M . 
More specifically, PSL(M|M) is obtained by first factor- 
ing out the U(l) subgroup thereby obtaining the sub- 
group SL(M\M) of GL(M|M), followed in a second step 
by the "gauging away" of the states capjdag the U(l) 
charges under j := J A A and j := J A j4 .r 3 tr 1 H 52 l This turns 
out to realize a line of RG fixed points (and confoemal 
field theories) labeled by the coupling constant g M E$ 



B. High-gradient operators when M, k —¥ oo 

We are going to show that, when k and M are very 
large, the spectrum for the one-loop scaling dimensions 



of high-gradient operators shares the same structure as 
that in Eq. (fh^). It will become clear by comparison to 
the case of k = 1 that the limit M, k — > oo is the extreme 
"classical" limit whereas the limit k = 1 is the extreme 
"quantum" limit. 

We restrict the family of "classical" high-gradient op- 
erators to objects of the form 



str (JJJJJJ) str (JJJ) 



(3.12) 



i.e., to diagonal GL(M|M)-invariant monomials of or- 
der s in both the holomorphic and antiholomorphic cur- 
rents. For any given order s, the engineering dimensions 
are all equal and given by 2s. This degeneracy is lifted 
to first order in the coupling constant g M . The task of 
enume ratin g all linearly-independent high-gradient oper- 
ators (3.12) of order s is greatly simplified by the as- 
sumption M, k — > oo. We can rule out the scenario by 
which it is a finite set of independent Casimir operators 
of g\(M\M) that fixes all the linearly independent clas- 
sical high-gradient operators of order s once the limit 
I\I — V oo has been taken. We can also rule out the sce- 
nario by which normal ordering changes the book-keeping 
between classical and quantum high-gradient operators of 
order s once the limit k — > oo has been taken. 

For high-gradient operators of type Eq. ( 3.1 0| ) or 
( |3.12 ) , the coupling g A does not renormalize their scaling 
dimensi ons , since g A (or O A ) can be removed from the 
action (|3.9[ ) by chiral transformation. All that therefore 
is needed to compute their one-loop scaling dimensions 
are their OPE with the quadratic Casimir operator O m . 
We will write the following expressions for the general 
case of GL(M|AT), and will set M = N (i.e., the case 
of interest) only in Eqs. (|3.19|), (|3.20j), and (|3.21|). The 



required OPEs follow from (a) the infra-trace formula 



str [JMJN] x str [Jj] 



= str ( JA/") str (M J) - str (M) str (jJAf) 

- str (JMJ) str (AT) + str (JM) str (JAf) 

(3.13a) 

and (b) the inter-trace formula 



str [jA^jstr [JJV] x str [jj] 



str (JJVJM) - str (JJNM) 
- str (JMAfJ) + str (JMJAf) 



(3.13b) 



with Ai and Af arbitrary operators. Here we have used 
the short-hand notation of Eq. (2.15). 

To proceed we also need to distinguish linearly inde- 
pendent high-gradient operators of order s. To this end, 
a "quantum number" , the number of switches, is intro- 



ducedfUij The number of switches of type and of 
type in a single trace are defined as follows. Consider 
the trace 



str (J„ 



1*2 Ma 



J, 



l>-2„ 



(3.14a) 
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where fi 1 , ■ ■ ■ , n 2n = ± while J_ = J and J + = J. Write 
the sequence of "conformal" indices 



Ml, 



! A t 2ni M2n+1 



(3.14b) 



where jjL 2r 



yLtj by cyclicity of the trace. The num- 



ber of switches of type j" is the number of sign 
changes from + — > — in two consecutive conformal in- 
dices when reading the sequence fi 1} ■ ■ ■ , /x 2n , M2n+i from 
left to right. The number n, of switches of type J, is 
the number of sign changes from — — > + in two con- 
secutive conformal indices when reading the sequence 

MlT ■ ' , M2n> M2n+1 Ir0m lcft to ri g nt - 

These quantum numbers are useful as it can be shown 
that there is no contribution in the one-lo op RG of super- 
traces made out of 2n currents as in Eq. ( 3.14a ) from the 



subspace with and to the one with at least + 1 
and n , + 1 . This implies a lower triangular structure for 
the linearized RG equati ons obeyed by all supertraces of 
order 2n as in Eq. ( 3.14a ) which allows to treat separately 
each sector defined by a given number of switches. We 
shall assume that the strongest renormalization of the 
engineering scaling dimensions occurs within the sector 
made of the maximum number of switches. 

Within the subspace of maximal switches it is sufficient 
to introduce 



u) := J, J — J J, 



fi m :=str( W m ), (3.15) 



is encoded by the operator 
R:=-2(N - M)J2 kn i 



(Z + n)0 ; O r 



In £1 



d d 



l+n 



l + n 



(3.18) 



It is instructiv e to compare the OPE (3.16) and the 
RG equation ( 0.18 ) with the corresponding result in 
the weakly coupled NLcrM on the symmetric space 
U(P + Q)/U(P) x U(Q) with P,Q > lffl They are 
essentially identical to the corresponding result for the 
U(P + Q)/U(P) x U(Q) NLcrM. 

Now we return to the case M = N. The diagonahaa- 
tion of R gives the largest and smallest eigenvaluesBEJJ 



Ail = +2s(s - 1) = -A« 



(3.19) 



Thus, both largest and smallest eigenvalues depend 
quadratically on s. In turn, one obtains a spectrum of 
one-loop scaling dimensions with the upper and lower 
bounds 



4L — 2s + s(s 1), 
9m 



» 



2s 



7T 



l -s(s-l), 



(3.20) 
(3.21) 



for any given 1 < s = 2, 3, • • • . Observe that these 
bounds are interchanged when g M — > —g M - 



for a ny to = 1, 2, 3, ■ ■ ■ . With the help of the OPE ( 0.1 3a| ) 
and (3.13b) one verifies the OPE 



stx{wu) m ujoj n ) x O 



M 



n+2 m 



str (ww m )str (ww") x O m 



-4n, 



(3.16a) 



m+n+2 ' 

str( WW m ) xO M = -fi^ — (N — M)O m+1 , 



and 



i 1 



fc+Z— m 

n m x O m = -2m ^ fi *A ~ 2m ( N ~ M )^r 

k,l=l 



n r ™ tt r n n xO M = -4r m r n mnn m+n n r ™ l Q r n ™ \ 

(3.16b) 

for any to, n, r m , r n = 1, 2, 3, • • ■ . 

The action of the linearized one-loop RG flow on the 
space of composite operators in the subspace of maximal 
switches spanned by 

L 

n?n?...fi?, $> r P = 2s ' ( 3 - 17 ) 

P =i 



C. High-gradient operators when k = 1 

Having dealt with the extreme "classical" limit, we 
turn our attention to the extreme "quantum" limit M = 
1,2, 3, • ■ • and k = 1 for which the interacting field the- 



ory (3.9) describes a problem of Anderson localization in 
d = 2 dimensions reviewed in Appendix [b|. 

The classification of all independent high-gradient op- 
erators in GL(M\M) or in PSL(M|M) is more involved 
than in SU(2) because the problem of listing all invariants 
is more complex.EiJ An increase of complexity can already 
be seen at the level of SU(iV) for which the invariant ten- 
sors of rank 2s are obtained from all possible products 
of one rank 2 tensor and two rank 3 tensors.Ej Instead of 
considering the most generic family of "classical" high- 
gradient operators (3.10), we consider the GL(M\M) in- 
variant family of "classical" objects 

{0S0£|m,n = O, 1,2,3,- •• , m + n = s} , (3.22) 

which must then be normal ordered. We are going 

( s) 

to prove that the coupling constant Zm'n of the high- 
gradient operator 0^0\ in the action 



2m V2tt 



m-\-n— s 



m.n— 



(3.23) 
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obeys the linearized one-loop RG equation 



j 7 0) 

dl 



(2-2s)Z^ n 



4^m(m - 1)ZW„ 

Z7T 



I — 5 

2tt 

4^(m+l)^« li „_ 1 



(3.24) 



for any m, n = 0, 1, 2, 3, ■ ■ • with m + n = s > 1. For the 
PSL(M\M) theory t he op erators 0^ are all absent. 

The RG equation (3.24) shows that there is no feed- 
back from high-gradient operators containing a factor 0\ 
to those containing a factor 0\ provided n' < n. Diag- 
onalization of the RG equation gives the set of one-loop 
scaling dimensions 



c W =2s + ^m( m -l) 



(3.25) 



for all m,n = 0, 1, 2, 3, • • ■ such that m + n = s. For a 
positive g M we get the lower and upper bounds 



41 = 2*, zLt = 2s + - 1), (3.26) 

7T 



respectively, i.e., Xm] n with 



n = s is always much 



larger than the engineering dimension 2 s so that the high- 
gradient operator O^O a is irrelevant. For a negative g M , 

(s) 

the spectrum of lower bounds on x m} n with m + n = s is 
unbounded from below when s — > oo, i.e., 



2s. 



(3.27) 



Proof: Having made the simplification g A = we 
only need to compute the OP E OV jO" x Q M where 
1 < m + n = s, t o jus tify Eqs. (|3.24|) and (|3.25|) . Each 
operator in Eq. (3.22) contains terms with 4s bosons, 
4s — 2 bosons and 2 different fermions, 4s — 4 bosons and 
4 different fermions, 4s — 2M bosons and 2M differ- 
ent fermions, and so on. The terms that contain identical 
fermions have short-distance singularities and hence they 
should be interpreted as operators that involve gradients 
over fermion fields after normal ordering. It is under- 
stood from now on that the OPE Oy^O\ x Om is only 
over the terms in the expansion 0^0\ involving differ- 
ent fermions, i.e., the OPE we present are "accurate" 
up to terms involving gradients over fermionic spinors. 
Neglecting the OPE between derivatives of the fermionic 
spinors and O m is harmless insofar as these OPE cannot 
feedback into the RG flows of those contributions that 
we keep. 

Let 



2M 



2M 



(3.28) 



A=l 



A=l 



and remember that Oa = — (V^W (V^Vy while Om = 
- (■0 t -0) {i>H)- The OPE that involve (ip^ip) and (ftip) 



(3.29) 



(^ t V>)(V> t V') x (l/)^1p)(^1p) 



-O 



A- 



On the other hand, the OPE that involve (tp^tp) , (V^VO > 
n/j''tp) ) and \M)'ip) are given by 



(^V)(^V) x (^VOW'V) 



(3.30a) 



M • 



(3.30b) 



Bot h C M and C A are generated through the OPE 
( 3.29) and ( 3.30 ), respectively. However, two OPE in 
Eq. (3.30a) always appear in a pairwise fashion and can- 
cel each other, 



(jj^){ fi))AxO M + (M)ty> T #4 x Q M = 0, (3.31) 

where „4 is some operator. Hence, the total number of 
O m contained in a high-gradient operator never increases 
under the linearized RG flow. 



From the OPE (3.2E) and (3.30) one deduces the OPE 



OmOI x O m 



mOZ^OlOu x Om + tiOZOI^Oa x O m 



+mnO M T 1 O n K - 1 OyiOA x O m 



+m(m - 1)0™- 2 OIOmO m X O m 



+n(n - \)O^QI- 2 Oa Oa x O m 
2m(m - 1)OmCa - 2m 2 0™- 1 O n + 1 



(3.32) 



(When m = 0, the term with 0^ is absent from the 



last line.) The linearized o ne-lo op RG equation (3.24) 
thus follows from the OPE (§]j§). □ 

Had we assum ed the level k to be larger than k = 1, 
the family ( 3.22j ) would not have been closed under the 



OPE with O m . For example, in the extreme classical 
limit M, k — > oo the family of hi gh-gr adient operators is 
given by the much larger family ( 3.12 ) . 

We close by pointing out that we could have reached 
the same conclusions on the spectrum of one-loop scal- 
ing dimensions of high-gradient operators had we used, 
instead of the effective action with diagonal GL(M|M) 
symmetry, an action built out of fermionic replicas or an 
action built out of bosonic replicas and taken the number 
of replicas to zero at the end of the day. Using bosonic 
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replicas mimics very closely the line of argument pre- 
sented here. Using fermionic replicas singles out high- 
gradient operators made of fermionic spinors that are all 
distinct through their replica index and then taking this 
replica index to zero, very much in the same way as repli- 
cated vortices in certain classeS|-af classical random two- 
dimensional Coulomb gasesE3~Ey 

We would like to stress that our results depend cru- 
cially on the continuous symmetry GL(2N\2N) of the 
g\(2N\2N) k=1 Thirring model. (From the point of 
view of Anderson localization, it is the existence of 
a continuous symmetry not the symmetry group per 
se that matters since the symmetry group changes de- 
pending on the choice made to represent single-particle 
Green's functions, say a supersymmetric, bosonic repli- 
cas, fermionic replicas, or Keldysh path integral.) If we 
consider local perturbations (local operators) that break 
the GL(2N\2N) symmetry, an infinite set of local op- 
erators with relevant (negative) scaling dimensions can 
appear. This alternative set of local operators may be re- 
lated to the-aituation recently considered by Le Doussal 
and SchehrEa The microscopic starting point of Rcf. ^ is 
a class of classical random XY models in two dimensions. 
These models can also be viewed as interacting models of 
Dirac fermions subjected to disorder, by the ma gi c o f the 
boson-fermion duality in d = (1 + 1) dimensions P 3 IP H The 
difference with our paper is that their model is not invari- 
ant under a continuous symmetry group, but only under 
the discretey symmetry group which permutes the replica 
indices. It is then necessary to use the full machinery of 
functional RG to account for the one-loop relevance of 
high-gradient operators. 



D. Comparison with the CP 1 ' 2 NLerM 

The perturbed gl(2N\2N) k=1 WZW model with g M > 
(Thirring model) describes a problem of Anderson lo- 
calization in two dimensions. As briefly reviewed in Ap- 
pendix this problem of Anderson localization arises as 
the long-wavelength description of a tight-binding model 
on a two-dimensional bipartite lattice with a form of dis- 
order that preserves sublattice and time-reversal symme- 
tries. The long- wavelength theory is a (2+l)-dimensional 
Dirac equation subject to disorder potentials consistent 
with these symmetries. In terms of the symmetry-based 
classification of Anderson localization, the relevant sym- 
metry jcjagis the class BDI (chiral-orthogonal symmetry 

It is possible to use a different representation of this 
Anderson localization problem, in terms of a NLerM 
whose target space is the non-compact supermanifold 



GL(2A|2A0/OSp(2iV|2A0. 



(3.33) 



model and the other in terms of the NLerM, are com- 
plementary to each other in that when one of the models 
is strongly coupled, the other is weakly coupled. A re- 
flection of this appears in the conductivity. The coupling 
constant of the NLerM is inversely proportional to the 
conductivity. In the clean limit g M = of the Thirring 
model the conductivity is of order unity (in units of e 2 /h) , 
consistent with the strongly coupled regime of the NLerM. 
The conductivity-increases with g M > as seen in per- 
turbation theory.Ej Furthermore, both g M and the con- 
ductivity are exactly marginal. This suggests a deeper 
relationship between the Thirring model and the NLerM, 
and indeed (following Ref. 41), one can turn the Thirring 
model into the NLerM continuously by tuning g M (or 
cquivalently the conductivity) continuously. 

We consider the case N = 1 for which the non-compact 



target supermanifold ( 3.33 ) is isomorphic to U(l) xU(l) x 
CP 1 ' 2 , where again a suitable analytical continuation is 
understood for CP 1 ' 2 , i.e., we need to consider the non- 
compact counterpart to CP 1 ' 2 as defined in Appendix |A|. 
Obtaining the non-compact CP 1 ' 2 target supermanifold 
of the NLerM from U(l) x U(l) x CP 1 ' 2 corresponds in 
the Thirring model to the reduction of the GL(2|2) to 
the PSL(2|2) current algebra in Eq. (|3 111). 

It is explicitly shown in Appendix |X| that all high- 
gradient operators are made more irrelevant at one-loop 
order by fluctuations in any non-compact QP N + M ^ 1 \ N 
NLerM labeled by the non-negative integers M and N. 
To be more precise, we find that the largest and smallest 
one-loop scaling dimensions for the high-gradient opera- 
tors of type ( |A14 ), for a given s, are 



.-<*> 



:2a + 2|i|a(s-l), 



a£L=2* + 2|t|xO, 



(3.34) 



(A suitable analytical continuation in the boson-bosen 
sector is needed to implement the non-compactness.E3) 
These two descriptions, one in terms of the Thirring 



where \t\ > is the coupling constant of the non-compact 
£pN-i\N NLc r]yL This is fully consistent with our finding 
( |3~26| ) in the Thir ring model. We conclude that, in sym- 
metry class BDI, high-gradient operators in the Thirring 
model with g M > behave in the same way as in the 
corresponding NLerM (i.e., the one that belongs to the 
symmetry class BDI). 

The sign of g M in the perturbed gl(2JV|2AT) fc=1 WZW 
model can be chosen to be negative, g M < 0. If so, 
this field theory does not represent anymore the mo- 
ments of the single-particle Green's function in a prob- 
lem of Anderson localization in (bulk) two dimensions. 
Nevertheless, this field theory does describe a problem 
of Anderson localization which, however, now belongs 
to the different symmetry class CII (chiral-symplectic 
symmetry class) describing the effect of disorder on the 
Dirac fermions which are known to appear at the two- 
dimensional boundary of a three-dimensionaLtopplogical 
band insu lator in the same symmetry class .e2tc3 Equa- 
tion (3.27) implies that high-gradient operators are now 
made more relevant by the current-current perturbation 
g M < to one-loop order. As for the case with g M > 0, a 
problem of Anderson localization in the symmetry class 
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CII is characterized by a NLctM with a corresponding tar- 
get manifold. As before, the beta function of the coupling 
constant g M < of the Thirring model as well as that of 
the coupling constant of the, corresponding NLcrM van- 
ish, and one can interpolatecil between the weak coupling 
limit of the Thirring model and the strong coupling limit 
of the NLcrM and conversely, by tuning g M continuously. 
The target supermanifold in symmetry class CII is the 
compact supermanifold (3.33), from which one extracts 
when AL = 1 the NLcrM on the compact supermanifold 
CP 1 ' 2 .^ It is explicitly shown in Appendix [A| that all 
high-gradient operators are made more relevant at one- 
loop order by fluctuations in any compact ^,p N + M - 1 \ N 
NLcrM labeled by the non-negative integers M and N. In 
particular, for M = 0, we find that the largest and small- 
est one-loop sca ling d imensions for the high-gradient op- 
erators of type (AAA), for a given s, are 



^ ax =2s-2ts(s-l), 



» 



:2s + 2t x 0, 



(3.35) 



where t > is the coupling constant of the com- 
pact CP^" 1 ^ NLcrM. Once again, we conclude that 
high-gradient operators behave in the same way in the 
Thirring model with g M < and in the corresponding 
NLcrM that belongs to the symmetry class CII. 



IV. CONCLUSIONS 

More than twenty years after their discovery, the role 
of high-gradient operators, which appear to be highly 
relevant in one-loop computations of anomalous dimen- 
sions in a great variety of NLcrMs, still remains a puz- 
zle. Indeed, this perturbative property is rather general 
as it can apply to both compact and non-compact tar- 
get manifolds. In the absence of an exact calculation of 
observables that would be sensitive to high-gradient op- 
erators, it is still an outstanding question whether the 
extreme RG-relevance of these operators is an artifact of 
the one- loop calculation (e.g., in the 2 + e-expansion) , or- 
is a feature that is generally valid. (For an attempt to 
compare the e expansion in d = 2 — e dimensions with 
exact results obtained for d = 1, see Ref. p^[ ) 

In order to shed some light on these issues we have 
asked in this paper the following question. Can high- 
gradient operators become relevant in the family of two- 
dimensional g\(M\M) k Thirring models with M and k 
positive integers? The strategy that we followed has 
three steps. The first step consists of identifying all the 
independent "classical" high-gradient operators of order 
s. This is a problem of group theory that involves the 
enumeration of all distinct GL(M|M) singlets in the di- 
rect product of 2s adjoint representations of GL(M|M). 
The second step consists of normal-ordering all indepen- 
dent classical high- gradient operators of order s. This 
step depends crucially on the level k of the non-Abelian 
Thirring model. The inverse level 1/k plays the role of 



a quantum parameter that vanishes in the limit k — > oo. 
The level k = 1 is thus the most "quantum". The com- 
putation of the linearized RG flows for the high-gradient 
operators is the final step. 

We could not solve the first step in its full general- 
ity. We were nevertheless able to construct two sets of 
high-gradient operators in the extreme "classical" limit 
g\(M\M) k with M, k — » oo and the extreme "quantum" 
case g\(M\M) k with M a positive integer and k = 1, 
respectively, and carry out the second and third steps 
consistently, i.e., show that each family of normal-ordered 
high-gradient operators is closed under the linearized RG 
flow equations. 

The set of high-gradient operators that we considered 
in the extreme "quantum" limit is much smaller than 
the set of high-gradient operators for the extreme "clas- 
sical" case. This is to be expected as normal ordering is 
extremely sensitive to the free-field fcrmionic representa- 
tion of the g\(M\M) k current algebra at the unperturbed 
WZW critical point. This difference has dramatic con- 
sequences for the spectrum of one-loop anomalous scal- 
ing dimensions in the extreme "classical" and "quantum" 
cases S3 

In the extreme "classical" case, anomalous one-loop 
scaling dimensions for high-gradient operators of order 
s are distributed in a symmetric fashion about zero 
with the minimum and the maximum both depending 
quadratically on the order s, very much like for the family 
of NLcrMs on the target spaces U(M-ftAUiU(A/) x U(7V) 
with M and N positive integers. E^Ej 123 Hence, high- 
gradient operators must become (one-loop) relevant for 
both signs of the current-current interaction with increas- 
ing order s very much in the same way as their cousins do 
in both the compact family U(M +N)/\5(M) x XJ(N) and 
the non-compact family U(M, N)/\J(M) x U(AT) with 
M, N > 1. 

In the extreme quantum case k = 1, the spectrum of 
anomalous one-loop scaling dimensions of order s is al- 
ways one-sided, i.e., positive for one sign of the current- 
current interaction. For g\(M\M) k=1 with M a posi- 
tive integer the sign of the current-current interaction 
for which high-gradient operators are always irrelevant 
corresponds to the interpretation of the gl(27V|2A r ) fe=1 
Thirring model as a problem of Anderson localization in 
random tight-binding models on two-dimensional bipar- 
tite lattices (symmetry class BDI) . We have shown in this 
paper that the high-gradient operators in these random 
tight-binding models are irrelevant at one-loop order. 

High-gradient operators in those NLcrMs of relevance 
to the physics of Anderson localization are related to the 
moments of the dc conductance.E3~E3 Their perturbative 
one-loop relevance has been interpreted as the signature 
of broad tails in the pro bability distribution of the con- 



ductance in Refs. 29 



(One should, however, bear in 
mind that the current-current correlation function enter- 
ing the Kubo formula for the conductance looks rather 
different from a simple GL(2A r |27V) current-current cor- 
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relation function. It would thus be very interesting 
to study the probability distribution of the dc conduc- 
tance in the relevant random tight-binding model using 
nonperturbative techniques (this may include, e.g., also 
numerical approaches) in order to establish if it is broad 
or not. 



Appendix [A| is organized as follows . W e first define 
NLcrMs on the projective supespaces ( |Al| ) using a geo- 
metrical approach. We then present the main result of 
Appendix A| on the one-loop scaling dimensions of high- 
gradient operators and show that they are one sided. We 
close by briefly outlining how the one-loop scaling dimen- 
sions of high-gradient operators are computed. 
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Appendix A: High-gradient operators in NLcrMs on 
the complex projective superspaces (£P N + RI ~ 1 \ N 

Whether or not the spectra of anomalous one-loop 
scaling dimensions of high-gradient operators in NLcrMs 
are symmetric about zero or not can be very impor- 
tant when the analytical continuation of the coupling 
constant t in the NLcrM from positive to negative val- 
ues is meaningful from a physical point of view. We 
shall call spectra which are fully symmetric about zero 
two-sided spectra. Spectra of anomalous one-loop scal- 
ing dimensions that are strictly positive (negative) will 
be called one-sided. NLcrMs with the target manifolds 



0(N)/0(N — 1) are already known to be one- 
sided.tJ We are going to show that thi&is also the case for 
the NLcrMs with the target manifoldsEO 



t£pN+M-l\N ^ 

\J{N + M\N)/[V(1) x U{N + M- 1\N)}. 



(Al) 



The complex projective superspaces (Al) are generaliza- 
tions of the compact complex projective spaces 



CP 



M-l 



U(M)/[U(1) x U(M - 1)]. (A2) 



We shall also study on their own right the high-gradient 
operators in NLcrMs with the non-compact complex pro- 
jective target spaces 

U(1,JV + M-1|#)/[U(1) x \J(N + M - 1\N)}. (A3) 

These non-compact manifo lds follow from the compact 
complex projective spaces ( |A2[ ) upon analytical contin- 
uation of some real coordinates to imaginary ones. The 
complex projective superspaces ( |A1[ ) are special cases of 
the Kahler supermanifolds whose high-gradient operators 
were studied in Ref. [l6| We refer the reader for notations, 
conventions, and the relevant intermediary results to Ref. 

m 



1. Geometry of the £pN+M-i]N NLcrMs 

A NLcrM on a Hermitian supermanifold can be defined 
with the help of the partition function 



Z := Jv[<p\<p] e- s ^>*\ 



(A4) 



Here, (ip^ip) are the coordinates on the Hermitian su- 
permanifold, a -gi,{Lp\ ip) is the metric on the Hermitian 
supermanifold, and t is the NLcrM's coupling constant. 

We are going to restrict our analysis to Hermitian su- 
permanifolds such that their metric can be derived from 
a Kahler potential. Furthermore, we shall choose the 
Kahler potential so that the target supermanifold is none 
but the CpN+M-iW symmetric space. The Kahler po- 
tential for CP N+M -^ N is 



K := log(l + <pi&). 



The bilinear form 



:= V a * ^ 



(A5a) 



(A5b) 



is presented in terms of the diagonal tensor £ with the 
components 

**ih~a*k (A5c) 
that do not depend on the coordinates (<p \ ip). Hence, 



1 + ^=1 + 



N+M-l N 



5>*V> (A5d) 



where (0**,^) with i = 1, ■ • ■ , N + M - 1 and 
with i = 1 , • • • ,N are the bosonic and fermionic co- 
ordinates of ((f'^ip), respectively. We observe that 
(tp* a , ip a ) has N+M—l bosonic and N fermionic coordi- 



nates. Equations (|Aj) and @ define the cpN+M-i\N 
NLcrMs. Setting N = in Eqs. (|A|) and ( |A5| ) defines 
the NLcrMs on the compact complex projective mani- 
fold (A2). The a nalyt ical c onti nuation cjf — > icj) 1 and 
</>** -> in Eqs. ([A§ and flAp ) defines the NLcrMs on 
the non-compact complex projective manifold ( |A3| ). 

The derivative of the Kahler potential ( A5a ) gives the 
metric (a superanalogue of the Fubini-Study metric) for 

C pN+M-l\N through 



~3 K d 

a-9b Qif a * dip b ~ a b 



(A6a) 
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Wc have introduced the scalar 



Z := 



1 



1 + f^(fi 
and the tensor 



(A6b) 



(A6c) 



Following the usage for graded indices from Ref. [l6|, we 
also define the tensors 



ayb* ._ a^b* +ip <^ (p b*^ 
\a-\-b-\-ab 



a Y b . ■■= (-1) 



,Y 



(A7) 



a*\/b 



ab b\/a* 



Y° - = (— i) a0 °Y 



a * 9b 


= z 




a9b* 


= z 


a^b* 7 


a *g b 


= z 


-la*yfe 




= z 





It then follows that the metric indices can be raised, low- 
ered, or shifted according to 



(A8) 



The metric tensor ( |A6| ) can be expanded about any 
point of the manifold at which it is finite (flat geometry). 
The lowest order in this expansion defines the "kinetic" 
contribution to the Lagrangian of the NLcrM, whereas 
the higher-order contributions define the "interactions" . 
The bosonic contribution to the kinetic ene rgy must be 
positive definite for the path integral (A4) to be well 
defined. This condition fixes the sign of the coupling 
constant t. For the compact complex projective mani- 



folds (A2), t > must be chosen. For the non-compact 
complex projective manifolds (|A3|), t < must be chosen 
(see, for example, Ref. |7(]). A consequence of the analyt- 
ical continuation i — >- — i for the one-loop beta function 
of t, if it is proportional to t 2 as it is in d = 2 dimen- 
sions, is that it changes by a sign. Similarly, the one-loop 
corrections to the scaling dimensions of high-gradient op- 
erators also change .by a sign under the analytical contin- 
uation t — ;> — t.E3E3 It then matters greatly whether the 
anomalous one-loop scaling dimensions are two-sided or 
one-sided. n 

According to Friedan,H the one-loop beta function for 
the coupling constant t of a NLcrM on a Ricmannian 
manifold is given by the curvature of the manifold. The 
curvature follows from the Ricci tensor, which we now 
compute for CP N+M -^ N . Needed arc the coefficients of 
the connection. They are 



pa — ad* 

1 be ■— 9 d*9b,c 



zK/V4 + (-i) bca <5 c /V4 



(A9a) 



and 



r a ; c . = a *g d d g b * 



(-l) ca 8 bV d d ^ + (-l) b+bca S c ip d d ^ ■ 

(A9b) 



The curvature tensor field on CP N+M 1 ' Ar can then be 
expressed solely in terms of the metric tensor field, 



bed* — 1 be° d 

= a 5 b ged*+ a 5e(-i) bc gbd 



and 



n b*c*d 



1 fc*c* ° d 



= a 5 b g c , d + a 8 c {-lf c g b . d . 



(AlOa) 



(AlOb) 



For CP N+AI 1 ' JV , the Ricci tensor field is proportional 
to the metric with M the proportionality constant, 



Rbd* = M gb 



(All) 



For (£p N + M ~i\N ^ f vj ows that the Ricci tensors van- 
ishes when M = 0, and so does the one-loop beta function 
according to Friedan. The beta function vanishes to all 
orders in the loop expansion. Ell The special case of CP 1 ' 2 
\(M, N) — (0, 2)] has also been discussed in Refs. ^ and 



High-gradient operators for the £pN+M-i\N 
NLcrMs 



From the property (A10), i.e., that the curvature 
tensor field of the supermanifold Cp Ar + A/ - 1 l w depends 
solely on its metric, follows that the RG equations among 
the infinite set of operators made of local polynomials in 



Guv : = d u ( P a *a'9b d v b( P 



(A12) 



are closed. 

Near two dimensions (d = 2 + e), it is convenient to 
use the conformal coordinates, 



d± = d x 



± id,, 



fJ, = ±,U = ±, 



(A13) 



i.e., we use G++, , Q h an d G-- as the building 

blocks for the high-gradient operators. It can be shown 
that the one-loop RG equations are closed within the 
family 



{Gl-G% (g ++ g..y }s:=i:= 



(A14) 



of high-gradient operators for any given number of gra- 
dients 2s, where p, q, and r are any non-negative integer 
satisfying p + q + 2r = s. Furthermore, for any given s, 



r, and r' the family (A14) obeys one-loop RG equations 
with an upper triangular structure in the sense that all 
high-gradient operators of the form (A14) with r' > r 
do not enter the one-loop RG equations for those high- 
gradient operators with r fixed. The task of diagonalizing 
the close d one-loop RG equations obeyed by the fam- 
ily (A14) thus simplifies greatly. It is indeed sufficient to 



fix s and r and to diagonalizc the one-loop RG equations 



obeyed by the family (A14) labeled by the non- negative 
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integers p and q. For any finite order s, diagonalization 
of the one- loop RG flows obeyed by the family ( A14 ) of 
high-gradient operators yields the one-loop RG eigenval- 
ues 



-2Mr 



+ 2(-pq + p(p-l) + q(q-lj), (A15a) 



here labeled by the non-negative integers q, q, and r that 
satisfy 



P- 



2r = s. 



(A15b) 



Combining Eq. ( A15 ) with the engineering scaling dimen- 
sion 2s yields the one-loop scaling dimensions 



x 



- 2s - tai s l 



0) 

p,q,r "^p^q^r 

2s - 2t 



^ — Mr —pq + P {p - 1) + q{q - 1)) 

(A16) 



for th e fam ily (|A14|) o f high-gradient operators. Equa- 
tions ( A15 ) and ( |A16| ) are the main result of this Ap- 
pendix. Observe that this result is independent of the 
integer N in CP iv + M ^ 1 l iv i Hence, it applies to the case 
iV = 0, b oth in its compact and non-compact incarna- 
tions (A2) and (A3), respectively. 

We now take a closer look at the spectrum when M = 
0. In this case, the projective supcrspacc is Ricci flat, i.e., 



R a , b = according to Eq. ( All ), and hence the one-loop 



beta function of the NLerM coupling constant t vanishes. 
(These areEJ in fact lines of critical points labeled by the 
coupling constant t of the CP w " 1|Ar NLcrMs.) We also 
distinguish the compact case from the non-compact case 
by demanding that t > in the former case and that 
t < in the latter case. 

The compact case corresponds to t > 0. For any given 
order s, we seek the largest and smallest one- loop RG 
eigenvalues that govern the RG flow of the high-gradient 
operators ( |A14[ ) in the NLcrMs CP™" 11 ™. Needed are 

(s) 

the extremal values of ctp, q , r while holding p + q + 2r = s 
fixed. We find that the most and least dominant one-loop 
scaling dimensions in two dimensions and for a fixed s are 



2s - 2ts(s - 1), 



£ S L = 2s - 2t x 0. 



(A17) 



We conclude that t he sp ectrum of one-loop anomalous 
scaling dimensions ( A15 ) for any "compact" <CP N ^^ N 
NLerM is one-sided in the sense that it is not symmetri- 

(s) 

cally distributed about zero: While x m [ n is not bounded 
as a function of s, XmLc = 2s irrespective of s. The re- 
sult for the most dominant scaling dimension x 
same as that for the 

U(P + Q)/U(P) x U(Q) 



0) 



is the 



(A18) 



NLcrMs with P,Q> However, the spectrum of one- 

loop anomalous scaling dimensions for the NLcrMs ( gig ) 



with P,Q > 1 is two-sided: The one-loop anomalous 
scaling dimensions are symmetrically distributed about 
zero. 

The non-compact case corresponds to t < 0. For any 
given order s, we seek the largest and smallest one- loop 
RG eigenvalues tha t gov ern the RG flow of the high- 
gradient o pera tors (A14) when M = 0. These follow 
from Eq. ( A17) with the substitution t —> —t, 



x£L=2s + 2\t\s(s-l), 
x ( ±=2s + 2\t\x0. 



(A19) 



So, there is no relevant high- gradient operator in this 
non-compact case. This is the conse quence of the one- 
sided property of the spectrum ( A15 ) when M = 0. On 
the other hand, in the case of the non-compact 



U(P,Q)/U(P) xU(Q) 



(A20) 



NLcrMs with P,Q > 1, there are always relevant high- 
gradient operato rs. We no te th at the one-loop scal- 
ing dimensions ( A17 ) and ( Al9j) turn in to th e corre- 
sponding scaling dimensions ( [3.27 ) and ( 3.26| ) for the 
gl(2AT|27V) fc=1 WZW model, if we identify t with -g M /ir. 

We now relax the condition for criticality M = 
of the CP N -^ N NLerM target mani fold. It can then 
also be shown that the spectra ( A15 ) labeled by s and 
M are one-sided. Since this result is, as required, in- 
dependent of N, it applies to the CP* 7 " 1 NLcrMs as 
well. In turn, CP M ~ 1 - U(M)/U(1) x U(M - 1) is 
obtained from U(P + Q)/U(P) x U(Q) by specializing 
to (P,Q) = (M - 1, 1) or (1, M - 1). The reason why 
the spectrum of one-loop anomalous scaling dimensions 
in U(P + Q)/U(P) x U(Q) with P, Q > 1 looks so dif- 
ferent from the cases with either P or Q being unity 
is the following. The U(P + Q)/U(P) x U(Q) NLcrMs 
with P, Q > 1 have a larger set of high-gradient oper- 
ators than in the projective (super) spaces. This can 
be se en b y comparing the set of high-gradient opera- 
tors ( A14 ) against their counterparts when the target 
manifold is U(P + Q)/U(P) x U(Q) with P,Q > 1, 
which can be found in Eqs. (2.12) and (2.16b) from Ref. 
£L§ High-gradient operators for U(P + Q)/U(P) x U(Q) 
NLcrMs with P, Q > 1 can be expressed product 
of traces of matrix fields, while in the complex projec- 
tive space, ther e is no such trace. [Here, note that Z 
defined in Eq. ( A6a ) is a scalar while the correspond- 
ing object in U(P + Q)/U(P) x U(Q) with P,Q > 1 is 
a matrix, Eq. (2.9b) from Ref. Similarly, the set of 
high-g radient operators in the gl fc>1 (M \M) WZW theory 
( |3.12 ) is larger than the set of high-gradient operators in 
the g\ k=1 (M\M) WZW theory. 



3. Sketch of the one-loop RG computation 

We now outline the calculations leading to the main 
results (Ejjjl) and flAl6|). 
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We choose the covariant background field method to 
renormalize the NLctMs. The merit of the background 
field method is that there is no need to worry about the 
appearance of redundant operators. This is very conve- 
nient when considering the mixing of a large set of oper- 
ators under the RG, that cannot be distinguished by the 
symmetries of the NLcrM. The background field method 
consists in resolving the coordinates (p a = (p^ + ( a of 
a NLcrM into slow (mean-field) modes ip^ that satisfy 
the classical equations of motion and fast (fluctuating) 
modes ( a in terms of which the Taylor expansion of the 
action transforms covariantly under reparametrization of 
the target manifold, i.e., in terms of which only the met- 
ric, the Riemann tensor, the covariant derivative, etc, of 
the target manifold appear in the action. For Kahler 
manifolds, this is achieved by choosing ( a to be (cither 
Ricmannian or Kahlcrian) normal coordinates. The very 
same expan sion of the action is also applied to the build- 
ing blocks ( A12) to the high-gradient operators, i.e., 



7 [IV 



[Guv] C o + [Q^v] 



c 1 



(A21) 



where [G^] ^ represents a p-th term in this expansion. 

To compute the anomalous scaling dimensions of high- 
gradient operators, they are first expanded in terms of the 
fast mode ( a , and are then pairwise Wick contracted. For 
example, the relevant formula for calculating ([(?^], 2 ) 

and ( [Gfj,u] ,i [Gpa] £i) s whereby the angular bracket (• ■ • ) 
denotes pairwise Wick contraction of the fast modes £ a , 



can be found in Ref. [T(| e.g., Eq. (C.40). When applied 
to the CP N+M - 1 \ N NLcrM, we obtain 

X [p^Gpu + G^Gpa) , 

(A22) 

where I = J d d k/{2ir) d (l/k 2 ). After substituting fi = 
±, v = ±, this gives 



<[0++] 


c-> = " 


-IMG+ 








([S-] 


c^> = - 


-IMG_ 








([G ++ ] 


c [ g +^ 




([G- + ] C [G + , 


W 


= 0, 


([S+-] 




-W = 


+4IG+.G+-, 






([3+-] 




-W = 


-2l(G ++ G^ 


+ G, 


-«-) 






-W = 


-2i(g__g ++ 


+ G- 


-+5+-) 


([G- + ] 













(A23) 

Furthermore, if p, q, and r are non-negative integers, wc 
find 



and 



{[Gl-G! 



>) = + P ([G + -] (2 )G p + Z 1 Gi + 
+ g^_<[S_ + ] c2 )^- + 1 
+ p<iGVgt + 1 ([G + J] c [G- + ] c ) 
+ P -^G^G1 + ([G + -] C1 [G + .] C1 ) 

+ q -^0l.Gl- + 2 ([G. + ] cl [G^ + ] cl ) 

= + 2l[ -pq + p(p - 1) + q(q - - 2Ipq G V +-G^G++G- 



(A24) 



[Gl_Gl + (G ++ G^)%) 



+ r(G p + .Gi + G r + - + 1 [G ++ ]^L_) 

+ r(G p + -Gl + G r ++ G-- 1 [V-] c >) 
+ ([Gl_Gl + ] c2 (G ++ G__) r ) 
-2IMrGl-Gl+(G ++ G—) r 

- 2l[pq-p(p - 1) - q(q - l)]G P + _Gi + {G ++ G^Y 

-2IpqG p + Z 1 G q r + 1 (G ++ G-j r+1 . 

I 



(A25) 



Equation (A25) justifies the claim that the family (A14) of high-gradient operators is closed under one-loop RG 
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and yields Eqs. (A15) and (A1G). 



Appendix B: Relationship between the gl(2A jr |2A) fe=1 
Thirring model and Anderson localization with 
"sublattice" symmetry 



a review 



The gl(27V|2A0 fe=1 Thirring model represents the 
physics of observables in a class of problems of Ander- 
son localization in symmetry classes BDI (see Ref. |33| ) 
and CII (see Ref. 41) within the classification scheme 
of Refs. [35| |38|. The fundamental physical observables 
are disorder averages of (products) of Green's functions. 
Here, we review some basic steps of this connection for 
the example of symmetry class BDI, whose simplest rep- 
resentative is a two-dimensional random tight-binding 
model for fermions on a bipartite lattice. A popular tight- 
binding model of that kind has recently become that on 
a honeycomb lattice, due to its relevance for the physics 
of graphene. (Another, lattice realization of the same 
continuum physics was obtained earlier in Refs. and 

Ho 



1. Definitions 

To be specific, consider the low-energy properties of the 
tight-binding model on the two-dimensional honeycomb 
lattice. Only sites on the different sublattices arc con- 
nected by hopping, and the hopping matrix elements are 
independent real random numbers with non-vanishing 
mean. Because of the constraint that only sites on differ- 
ent sublattices are connected, the model inherits a spe- 
cial symmetry called sublattice (or chiral) symmetry. It 
turns out to imply the presence of an operator that anti- 
commutes with the Hamiltonian, which thus relates the 
spectrum at positive and negative energies, and makes 
the zero of energy, E = (often called the "band cen- 
ter" ) , special. Taking the low-energy limit near zero en- 
ergy one obtains a random Dirac equation. (See, e.g., 
Ref. |5| for details.) We now start from the continuum 
limit of the so-obtained Hamiltonian, which reads 



H = H Q + V(r), 



(Bla) 



where the kinetic energy is (we set H and the Fermi ve- 
locity v F to be one) 

2 

^ = -E z K® r i)^> ( Blb ) 

and the static disorder is 
2 

V(r) = ^K® T 2)^( r )-( cr o® r 2 )V{r) 
+ {<73®n)M(r). 



Here, <r 1 2 3 and t x 2 3 are two independent sets of Pauli 
matrices together with another two independent 2x2 
identity matrix er 



I 2 and r 



I,. The 2x2 ma- 



trix space associated with the r Pauli matrices origi- 
nates from the bipartite symmetry of the underlying lat- 
tice model. The real- valued functions (potentials) A^r), 
V(r), and M(r), which do not vary appreciably on the 
scale of the lattice spacing, represent four independent 
sources of (static) randomness. 

The above potentials are random variables. We will 
assume first that they are white-noise distributed accord- 
ing to a Gaussian probability distribution with vanishing 
mean. 



A^(r)AJr>) = g A 5^ 2 \r - r% fx,u=l, 2, 



V{r)V{r') = M(r)M{r') = g M 5 {2) (r - r'), 



(B2) 



where 5^{r — r') is the two-dimensional delta function, 
777 represents disorder averaging, and we assume that the 
variances of V(r) and M(r) are identical. Of course, the 
disorder strengths g A and <? M are positive. 

The tight-binding Hamiltonian is invariant under time- 
reversal and so is its continuum limit 



T(U)*T = H, T: 



(B3) 



Since the tight-binding Hamiltonian preserves the bipar- 
tite nature of the underlying lattice for any realization 
of the disorder, so does its continuum limit through the 
chiral symmetry 



C1~LC — — 7^, 



C := <7n 



(B4) 



As already mentioned, because of its chiral and time 
reversal symmetries the Hamiltonian belongs to the 
so-called BDI s ym metry class within the classification 
scheme of Refs. B5H38. 



2. Path integral representation of the 
single-particle Green's function 

In problems of Anderson localization, physical quanti- 
ties are expressed by disorder averages of (products of) 
the retarded and advanced Green's functions 



G R/A (E) := (E±ir 1 -H)- 



(B5) 



In the present model, the retarded and advanced Green's 
functions are related, at the band center E = 0, by the 
chiral symmetry through 



C G R (E = 0) C = —G A (E = 0). 



(B6) 



Hence, any arbitrary product of retarded or advanced 
Green's function at the band center equates, up to a sign, 
a product of retarded Green's functions at the band cen- 
ter. From now on we will omit the energy argument of 
the Green's function having in mind that it is always 
fixed to the band center E = 0. 
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Since the two kinds of Green's functions are related by 
the chiral symmetry, it suffices to introduce functional 
integrals for the retarded Green's function defined by the 
supersymmetric partition function Z = Z F x Z B withu 



Z F := J V[x, x] exp \ i J X (iV - X 
Z B ■= I V[U]e*p(i [ £(ir)-H)Z 



(B7) 



Here, J r = J d 2 r = J dzdz/(2i), (x, x) is a P a i r °f two in- 
dependent four-component fermionic fields, and (£, £) is a 
pair of four-component bosonic fields related by complex 
conjugation. 

The matrix elements of the retarded Green's function 
can be represented as 



iG R (r 7 r')=( X (r)x(r')) 
=«(r)f(rO> 



(B8) 



with (• • • ) denoting the expectation value taken with the 
partition function Z. With the help of the property T — 
T 1 ' of time-reversal and the property in Eq. (B3), 



f X(ir) -H)x = ~ [ X T T{ii! - H) T X T , 

J r J t 

aiv - u)t = + [ er(t V - u) tf, 



(B9) 



Eq. ( p38|) is also given by 



iG R (r,r>)=-((T X T )(r)(x T T)(r')) 

= + ((rc T )(r)(e T r)(r')). 



(BIO) 



We now perform the change of integration variables 
Xi X ~^ V* ) "0 an d — ^ P iP where, when the matrix 
space on which the r Pauli matrices acts is made explicit, 



\l ^{^Vxi-i^x) i X 



I 

■la 
1 



+# 2t 
Pi 



We also define 

A := A„ + iA 



yl ■ "i yi 

m :=V — iM, m := V + iM, 



(Bff) 



(B12) 



With these changes of variables, the partition function at 
E = can be written as 



-S* -Sf. 



Z B = I V\p\p] e- sB - s > 



(B13a) 



with the effective action for the fermionic part given by 
[V> at (23 + A)^ a + ^ (2d + A)^ a 



2tt ^ 

a=l 



"t, 



(B13b) 



q F 

ir] 



and the bosonic part of the effective action given by 



(B13c) 



Observe the non-Hermitian appearance and asymmetry 
between fermions and bosons in j2fZ5 , which are neces- 



sary to maintain supersymmctryt 



The time- reversal invariance (B9) and 
of the new basis implies invariance underL. 



fO) in terms 



V>2 



V>1, 



-iPi, 



Pi ->- 



(B14) 



The finite level-broadening term ' is necessary for 
the computation of certain physical observables, includ- 
ing for example the Kubo conductivity, the Einstein con- 
ductivity, and the local density of states. O However, 
when we compute the conductance from the Landaucr 
formula by attaching ideal leads to t he disordered region 
described by the Hamiltonian ( Bf a ), we can set r\ = 
in the .disordered region (while still keeping r\ ^ in the 
leads) M 

The last step consists in averaging the partition func- 
tion Z = Z F x Z B over the probability distribution for 
the white-noise and Gaussian distributed random poten- 
tials. In this way, one finds a generating function for the 
averages of Green's functions which is nothing but the 
gl(2[2) fc=1 Thirring model. Specially, integration over the 
vector potential yields the term proportional to 



=(-) X ^ At x (~) B M m , (B15a) 



while integration over the complex- valued mass yields the 
term proportional to 



4> B H B x i> A ^ A =(-) A ^ A ^ ^ B ^ A \ 



(B15b) 



where we have combined bosonic /3 a , /3 a and fermionic 
i> a ,ipa (a = f , 2) spinors into the supersymmetric vector 
ip A , ip A (A = 1, • • • , 4) as in Sec. HI. The iV-th moment 
of the retarded Green's function evaluated at the band 
center is obtained by allowing the index a to run from 
1 to 2N in Eq. (BI3) or, equivalently, by allowi ng th e 
indices A and B to run from 1 to 2N + 2N in Eq. flBI5| ), 
thereby obtaining the gl(2N\2N) k=1 Thirring model. 
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We have assumed so far that the random imaginary cumulants (but still assuming that they have no spatial 

vector potential A, A and the complex random mass m, m correlations), the quenched disorder averaging necessar- 

possess a Gaussian probability distribution. If we assume ily yields high-gradient operators of the form ( 3.22 ). 
instead that their distributions have non- vanishing higher 
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